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Abstract. 

Spectropolarimetry is a photon-hungry technique that will reach fruition in the 8-m telescope age. 
Here I summarize some of the stellar spectropolarimetric research that my collaborators and I have 
undertaken, with particular emphasis on the circumstellar environment of massive stars, symbiotic 
binaries, and star formation. 



INTRODUCTION 

Stellar spectropolarimetry 1 is a relatively underused technique that has enormous di- 
agnostic potential, particularly for research into circumstellar matter and close binary 
systems. The principle difficulty with the method is that very high signal-to-noise ob- 
servations are required, since typical astrophysical polarimetric signatures may be on 
the order of a few tenths of a percent. When photon noise is dominant, the polarimetric 
precision AP goes as 

100 

AP^— (1) 



N 

where AP is in percent and N is the number of photons. One can see that a S/N of 1000 
or more per pixel is needed to measure a spectropolarimetric signal at the 0.1% level. At 
moderate resolutions (R ~ 5000) this limits even 4-m class telescopes to relatively bright 
objects (say V < 12). Higher dispersions and fainter objects necessitate even larger 
apertures, and fortunately spectropolarimetric instrumentation is present or planned on 
most of the world's 8-m facilities. 

The origin of linear polarization in a stellar context is usually dichroic absorption 
by dust, or scattering by dust, molecules or electrons (although linearly polarized ther- 
mal emission from aligned grains is sometimes observed in the millimetre regime). For 
an unresolved source an asymmetric geometry with respect to the observer's sightline 
is required, for example the alignment of grains by a magnetic field, an asymmetric cir- 
cumstellar distribution of material such as a disc, or intra-binary scattering. Cancellation 
between polarized radiation from different regions often reduces the observed polariza- 
tion to the level of a few percent - nonetheless the polarimetry enables one to probe 
geometrical parameters of an unresolved object. 



Note that this paper focuses on the use of linear spectropolarimetry, although of course circular spec- 
tropolarimetry is an extremely powerful tool for probing stellar magnetic field structures, both at the stellar 
surface and in circumstellar material (e.g. GlBl)- 
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FIGURE 1. Spectropolarimetry of the 6825A Raman-scattered line of D-type symbiotic RR Tel. This 
'triplot' shows the intensity spectrum (bottom panel), the polarization magnitude in percent (middle 
panel), and the position angle of the polarization (top panel). 

In this paper I present a synopsis of some of the stellar spectropolarimetric studies 
that my collaborators and I have undertaken. I concentrate on three broad areas: Raman- 
scattering in symbiotic binaries, and its use in measuring orbits in these long-period 
systems; the use of recombination-line polarimetry to measure global structure in the 
winds of massive stars; and the use of spectropolarimetric in examining the circumstellar 
structure of pre-main- sequence stars. In the final section I give a brief account of the 
instruments and likely avenues of research in the 8-m telescope age. 



RAMAN SCATTERING IN SYMBIOTIC BINARIES 

The spectra of symbiotic stars exhibit both molecular absorption bands and nebular 
emission lines. It is now known that this so-called 'combination spectrum' is formed due 
to binary interaction between a mass-losing red-giant star and a hot companion (typically 
a white dwarf), which ionizes the cool wind from the giant. Symbiotic stars are thought 
to be detached systems Jll], with accretion onto the hot component occurring from the 
wind rather than via Roche-lobe overflow. Systems that show a stellar spectrum in the 
IR are classified as S-type systems, while those that show dust emission are denoted 
D-type. 

In a seminal paper Schmid |4] determined that the AA6825,7082A emission lines 



0'75. 




1 0.5 -0.5 -1 

FIGURE 2. The polarimetric orbit of SY Mus. The O VI source is located at (0,0) and the filled dots 
show the position of the red giant every O.Olf. Short thick 'bowties' indicate the measured PAs and 
their uncertainties, the long thin lines show the calculated PAs at the times of the observations. The 
spectroscopic phases of quadrature and conjunction are also marked. 



that are observed in ~50% of symbiotic systems arise from Raman scattering of the 
A A 1032, 1038A O VI resonance doublet by neutral hydrogen. Raman scattering is the 
inelastic analogue of Rayleigh scattering: the frequency of the scattered photon is dif- 
ferent to that of the incident photon, and the scatterer is left in an altered quantum me- 
chanical state. In this case the O VI photons scatter in the far-red wing of Ly/3 and the 
scattering H° atom is left in the n = 2 state. By conservation of energy the frequency of 
the Raman photon, v r , is given by 

Vr = V P - V Lya (2) 

where VL y « is the frequency of the Lya transition and Vp is the frequency of the parent 
photon. 

The asymmetry of the scattering geometry, with O VI photons produced near the hot 
star scattering in the cool star's wind, coupled with the dipole nature of the scattering 
event, leads naturally to linear polarization [5]. Extensive surveys of symbiotic systems 
(d&LZD) revealed that the Raman lines often have a complex structure, typically compris- 
ing two or three peaks in both intensity and polarization accompanied by a 90° position 
angle (PA) flip and/or PA rotation (see, e.g., Figured]). The symbiotics show large line 
polarizations (up to ~ 15%), with the strongest and most highly polarized lines found in 
D-type (dusty) Mira systems. 

Numerical modelling of the Raman-line formation process (10, EE Ell) indicates that 
the PA of the blue wing of the Raman-line polarization is perpendicular to the plane 
containing the binary line-of-centres and the observer's line of sight. The models further 
demonstrate that the magnitude of the line polarization (P/) is also a simple function of 
orbital phase 

P,=Asin 2 a (3) 
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FIGURE 3. Plots of (a) orbital period against spectral type, and (b) the distance to the L\ point against 
stellar radius. The targets from this study are plotted as filled squares, and the open circles correspond 
to other symbiotics of 'known' period listed in |3]. The dashed line in (a) corresponds to the minimum 
period relationship while the dotted line connects the longest period systems. The dotted lines in panel (b) 
correspond to R = L\ , R = 2L\ , and R = 3L\ . The solid line joins the maximum radius of CD— 43° 14304 
(1OOR ) to the radius predicted from its spectral type (48R ). 



where a is the angle between the observer's sightline and the binary line-of-centres, and 
A is a scale factor which accounts for geometrical and physical parameters (interbinary 
distance, stellar radii; mass-loss rate, velocity law, etc). By observing the change in line 
polarization with time, it is possible to derive orbital parameters (e.g., period, inclination, 
PA of line-of-nodes) with analogous methods to those used for visual binaries. The 
first 'spectropolarimetric orbit' ,that of SY Mus (Figure El 111 IIP , was followed by 
those of AG Dra and Z And (El [l3]), while follow-up observations yielded orbits 
of four additional systems (CD-43° 14304, Hen 1242, Ml-21, V455 Sco) JH, thereby 
doubling the number of known orbits with P > lOOOd. 

It has been shown by Miirset & Schmid ([3]) that there is a strong correlation between 
orbital period and spectral type of the cool component, with the majority of systems 
lying close to a straight line relating the spectral subtype to the orbital period. For a 
given subtype no systems are found with substantially shorter periods. The minimum 
period locus is given by 

log/ 5 min = 0.1175 + 2.28 (4) 

where S is the M-class subtype and P m [ n is in days (see Figure |3h). Miirset & Schmid 
proposed that this period limit is related to the red giant radius, which was found to be 



L\/2 for most systems (where L\ is the distance from the primary to the inner Lagrangian 
point). Essentially this means that for binary to appear as a symbiotic it must be close 
enough for significant wind accretion onto the hot component, but wide enough that the 
binary is detached (Figure |3f>). Our systems with spectropolarimetric orbits conform to 
this relation, excepting CD— 43, which has an exceptionally early spectral-type for its 
cool component (K7). 

The spectropolarimetric approach to orbit determination has unique advantages. First, 
because the observed quantities can be determined to a small fraction of the orbital 
amplitudes (e.g., PA to ~ 1° in 180°), the approach is extremely efficient - e.g., in 
the favourable case of SY Mus, we were able to determine correct orbital characteris- 
tics from just four observations, whereas the same information from traditional radial- 
velocity techniques required many observing nights |fl5ll . Secondly, the amplitude of 
variation of the polarization is independent of orbital period, so that adequately sampled 
datasets are capable of yielding accurate orbital characteristics for systems with periods 
of from a few months to many years (unlike normal spectroscopic methods). Our group, 
along with Schmid's at ETH Zurich, are continuing to obtain data on various symbiotics 
with the aim of establishing periods for many more S-type systems, and obtaining the 
first reliable D-type orbit. The latter goal is obviously long-term, since D-type systems 
are likely to have periods of many decades ! 

THE GEOMETRY OF STELLAR WINDS 

Massive stars shed mass via high-speed, radiation-driven winds. Mass-loss rates of 10 -7 
to 10~ 4 MQyr" 1 are typical, with terminal speeds of 1000-3000 km s -1 . Stellar wind 
models are often used to fit recombination line profiles (such as Ha) in order to deter- 
mine mass-loss rates. This approach essentially reduces to 'counting' recombinations, 
based on the standard assumptions of homogeneity and sphericity. Departures from ei- 
ther will lead to an overestimate of the true mass-loss rate: a wind that is clumpy on 
small scales, or is global distorted, will have an increased rate of recombinations (which 
is a density-squared process) per unit mass, so a structured wind will produce greater 
line emission than a smooth wind with the same mass-loss rate. 

So how applicable is the assumption of sphericity? The winds cannot be imaged 
directly, but fortunately polarimetry provides an indirect probe of the geometry of the 
wind near to the star. Continuum photons from the star will pass through the ionized 
wind, and some will be electron- scattered which introduces linear polarization to the 
starlight. The nett polarization will be zero for a spherical distribution of electrons, but 
will be non-zero for any geometry that is not circularly symmetric with respect to the 
observer's line of sight. The recombination line emission, which is produced further 
out in the wind, will 'see' a lower electron-scattering column, and will therefore be 
less polarized than the continuum. Spectropolarimetry of a non- spherical wind should 
therefore show a polarized continuum (on the order of a few tenths of a percent) with the 
magnitude of the polarization dropping at the emission line wavelengths, as the polarized 
continuum is diluted by unpolarized line emission. This known is the 'line effect' . As an 
illustration I have plotted polarization spectra of two Wolf-Rayet (W-R) stars in Figure|H 
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FIGURE 4. Spectropolarimetry of the WN star WR 134 and the WC star WR 137, binned to a 
constant polarization of 0.1% and 0.05% respectively. Note the strong depolarizations at the emission 
line wavelengths. 



note the obvious decline in continuum polarization magnitude across the emission lines, 
a sure indication that the stellar winds in these objects depart from spherical symmetry. 

For WR137 we have other indicators that wind may be asymmetric. This WC star 
has an O-star binary companion in an eccentric orbit with a period of ~ 12 years. At 
periastron the interaction of the winds from the two components provides the necessary 
conditions for the condensation of warm carbonaceous dust grains, which produce a 
near-IR excess iflfjll . There is an intimate relationship between the solar-radius-scale 
geometry of the stellar wind from polarization, and the larger scale distribution of newly 
formed dust (many hundreds of AU) - see Figure |5] and [17]. The link between near- 
star geometry and a large-scale nebula has also been established for the Luminous Blue 
Variable (LBV) AG Car lfl8ll . Spectropolarimetry of this object shows a strong line effect 
at Ha, along with substantial variation (> 1%) in continuum polarization. The PA of the 
intrinsic polarization is aligned with the major axis of the AG Car ring nebulae and 
perpendicular to its jet, indicating that the resolved circumstellar environment is already 
present within a few stellar radii of the central object. 

The magnitude of the continuum polarization (P c ) for a given stellar wind geometry 
and viewing angle ( i) may, in the optically thin limit, be computed analytically 12011 : 

P c «2r(l-3y) sin 2 / (5) 

where T is the angle- averaged electron-scattering optical depth of the envelope and y 
is a shape factor. The degeneracy between the stellar wind geometry and the viewing 
angle means that studying objects on an individual basis is rather difficult, but if a large 
enough sample is obtained it is possible to disentangle the two statistically. We gathered 
a spectropolarimetric data set of 16 W-R stars using the 4-m William Herschel Telescope 
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FIGURE 5. The 1998 K' maximum entropy restored images of WR137 (logarithmic greyscale), and 
the 1997 image (contours) adapted from Fig. 3 of 1 19]. The dashed line illustrates the PA of the intrinsic 
polarization. 

(WHT) and the ISIS spectrograph JH. A statistical analysis of these data using Monte- 
Carlo based K-S tests enabled us to show that about 15% of W-R stars show a global 
departure from spherical symmetry. The level of the intrinsic polarization is consistent 
with an equator-to-pole density ratio of 2-3. We postulated that the line-effect systems 
represented the most rapidly rotating W-R stars, in which the rotation was having a 
significant affect on the wind dynamics. 

Interestingly the line-effect systems are clustered towards high mass-loss rates and 
high-luminosity (Figure |6l). Although at first sight this appears to suggest that this may 
be the root cause for the wind asymmetry, one should note that these quantities are very 
dependent on small and large-scale structure, as noted above. In other words the high 
mass-loss rates and luminosities determined by the models may be a result of applying 
spherically symmetric models. Mass-loss rates derived from thermal radio emission, 
which are the least model-dependent, show good agreement with those derived from the 
optical recombination lines. This does not mean that the mass-loss rate values are correct 
in an absolute sense, since both the radio continuum and the recombination lines scale 
as density-squared, but the agreement does indicate that the degree of wind asymmetry 
is similar at the radii for optical depth ~ 1 of both the optical ~ lO 1 ^* and the radio 
~ 10 3 tf*. 

Since departures from spherical symmetry appear to be relatively common for the W- 
R stars, and also occur in LBVs, it seemed natural to search for evidence for global struc- 
ture in the winds of O-supergiants. We obtained Ha spectropolarimetry of a sample of 
20 Northern, and combined these data with AT-band spectropolarimetry obtained with the 
4-m United Kingdom Infrared Telescopes and low resolution optical spectropolarime- 
try from Pine Bluff Observatory lE2fl . Somewhat surprisingly, the spectropolarimetry 
revealed that the single O-supergiants typically have intrinsic polarizations of < 0.1% 
(which corresponds to an equator:pole asymmetry of approximately less than 1.25). It 
appears that the wind line-driving mechanism in rapidly-rotating O supergiants (whose 
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FIGURE 6. Wolf-Rayet stars plotted in the mass-loss/luminosity plane. WN stars are plotted as squares, 
while WC stars are plotted as triangles. Stars with a line effect are displayed as open symbols, stars with 
no detectable line effect are shown as filled symbols. 

winds are optically thin) is relatively insensitive to rotation compared to the optically- 
thick line-driving mechanism of W-R outflows. 

W-R stars are also known to have small-scale clumps that appear to accelerate slowly 
from the star (e.g. ||23|]), and there is spectroscopic evidence for similar structures in 
the O supergiant £ Puppis [24]. These migrating density structures could potentially be 
identified polarimetrically, but numerical models of density enhancements propagating 
through the wind predict polarizations of close to 0.1%-0.2% [25], which is uncomfort- 
ably close to the current limits on the absolute precision of polarimetric observations. It 
is clear that very precise techniques will be needed to probe these clumps in the winds 
of O supergiants. 

PRE-MAIN-SEQUENCE STARS 

The broad picture of low-mass star formation, from the collapse of a gravitationally 
unstable molecular cloud, through to a pre-main-sequence star surrounded by a dusty 
accretion disc, is well-established. During the final, classical T Tauri star (CTTS) phase 
magneto spheric accretion is thought to occur, with the star's magnetic field disrupting 
the disc and channelling the material onto the surface, where its kinetic energy is liber- 
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FIGURE 7. Ha spectropolarimetry of a Herbig Be star (HD 53367) and a Herbig Ae (XY Per), binned 
to a constant error in polarization of 0.03% and 0.1% respectively. Note that the Herbig Be star shows a 
classical line depolarization, whereas the polarization spectrum of the Herbig Ae star is more complex. 



ated at hot spots. However, the situation for intermediate and high-mass star formation is 
much more confused, with a lack of clear-cut evidence for circumstellar discs. In the pre- 
vious sections I have shown that line spectropolarimetry can be a powerful probe of the 
circumstellar geometry, and we are now in the process of applying spectropolarimetric 
techniques to pre-main-sequence (PMS) stars. 

We have obtained a large sample of Ha polarization spectra Herbig Ae/Be stars lEfill . 
The two classes of star show a dichotomy in spectropolarimetric characteristics. The 
Herbig Be stars typically show the classical depolarization line effect, which is indicative 
of a large volume of Ha emitting material embedded in an asymmetric envelope. The 
Herbig Ae stars on the other hand show a more complex polarization through the line, 
with some stars showing significant polarized Ha emission. This suggests that Herbig 
Ae stars have compact sources of Ha emission, possibly magneto spheric accretion 
streams, with Ha photons scattering in the surrounding asymmetric medium. It is not 
completely clear whether the scatterers here are dust or electrons, although thermal 
motions of electrons should lead to significant smearing of the line polarimetry, which 
is not observed. Our results indicate that the magneto spheric accretion paradigm for 
low-mass star formation may extend to higher masses. 

A significant investment has been made in photopolarimetric studies of classical 
T Tauri stars (CTTS), which may show continuum polarizations of a few percent e.g. 
|27l|28l|29l|30]. This polarization is variable, perhaps due to changes in the illumination 
of the scattering envelope due to hot spots 113 ill . Since it was thought that the line and 
continuum photons from the central object were scattered by the same dusty envelope 
the expectation was that there should be no line effect. However our recent observations 
of RY Tau (Figure HI [32]) show a dramatic change in the polarization through Ha. It 
appears that RY Tau is not exceptional in this regard: in our recent survey 9 out of 10 
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FIGURE 8. Ha spectropolarimetry of the classical T Tauri star RY Tau, binned to a constant error in 
polarization of 0.1%. 



CTTS showed a detectable change in polarization across Ha Q. These observations 
are qualitatively consistent with scattering of a compact source of Ha photons off a 
circumstellar disc. 

A more detailed interpretation of the spectropolarimetry requires numerical mod- 
elling. As a first attempt we have examined the polarimetric signature of dipolar scatter- 
ing of a compact line source off optically-thin and optically-thick dusty discs [34] using 
the Monte-Carlo radiative-transfer code TORUS W25I1 . We find that the PA dependence 
of the line profile is an indicator on the extent of the disc, with a double PA rotation 
through the line corresponding to a disc with a negligible inner hole, and a single PA 
rotation indicating a substantial evacuated area around the central star. By interpreting 
the observations in the context of these numerical simulations we will begin to place 
strong constraints on the circumstellar environment of PMS across the mass spectrum. 



THE FUTURE 



Spectropolarimetric capabilities are already available on 8-m class telescopes. For ex- 
ample, VLT/FORS offers both circular and linear spectropolarimetry, but at R ~ 2500 
its maximum resolution is insufficient for most of the line polarization work described 



here (although the broad emission lines in supernovae can be resolved). Similarly the 
Keck/LRIS combination is limited essentially to continuum studies and is used predom- 
inantly for extragalactic research. The two 8-m Gemini telescopes are due to be fitted 
with facility polarization modulator units (GPOL), which will eventually provide a num- 
ber of their optical and near-IR instruments, including a high-dispersion echelle spectro- 
graph, with a polarimetric capability. The Prime Focus Imaging Spectrograph (PFIS) on 
the 10-m South African Large Telescope will provide a medium resolution (R ~ 5000) 
spectropolarimetric capability, and is ideally suited to stellar work. 

Pleasingly it appears that spectropolarimetry is now viewed as a standard mode by 
spectrograph developers. The proliferation of spectropolarimetry in the 8-m telescope 
age will lead to exciting new advances in stellar astrophysics, for example in the in- 
vestigation of the link between rotation and outflow asymmetries and structure in mas- 
sive stars and asymmetries in core collapse supernovae, which can show a line effect 
( l35ll36ll ). and may have viewing-angle-dependent observed luminosities. Rotation and 
asymmetry obviously play a key role in the collapsar model for gamma ray burst progen- 
itors tyi\\ . and high-resolution linear spectropolarimetry of the burst afterglows should 
yield valuable insights. Extending the spectropolarimetric study of W-R stars into the 
low-metallicity environment of the Magellanic clouds would provide a good test of line- 
driving theories. The move to longer wavelengths will also lead to new discoveries, for 
example by probing the circumstellar environments of highly obscured objects, such as 
Class I protostars that are still embedded in their nascent molecular clouds. 
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